Global singularity-free solution of the Iordanskii force problem 
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We present a derivation of the transverse force acting on a hydrodynamic vortex in the presence 
of an incoming sound wave from a global solution of the scattering problem, using the method 
of matched asymptotic expansions. The solution presented includes a detailed treatment of the 
interaction of the incident wave with the vortex core, and is free from the singularities in the 
momentum exchange between vortex and sound wave which have led to contradictory results for 
the value of the transverse force in the literature. 



I. INTRODUCTION 

The transverse force on a quantized vortex due to 
phonon scattering, first found by Iordanskii Q, has been 
controversially discussed in particular since the work of 
Thouless, Ao and Niu, which predicted that it in fact 
would vanish @, USEES HUGHEO] Based 
on general properties of superfluid order and the Berry 
phase associated with the vortex motion, they came to 
the conclusion that the total transverse force is propor- 
tional to the superfluid density alone and given solely by 
the superfluid Magnus force, and not, as the existence 
of the Iordanskii force would imply, proportional to the 
total density, normal plus superfluid. Recent calculations 
based on the vector potential due to the Berry phase of 
quasiparticles lead to the same conclusion of a vanishing 
Iordanskii force (Tll |. 

It therefore seems in order to reinvestigate the as- 
sumptions underlying calculations of the Iordanskii force 
based on classical (but two-fluid) hydrodynamics (see, 
e -g-> PL E|); where the transverse force stems from the 
momentum transfer between asymmetrically scattered 
sound waves and the vortex [l2L Il3l | . The existing treat- 
ments use various simplifications for a treatment of this 
problem, which constitute sources for possible ambigui- 
ties and problems in the final result obtained for the value 
of the transverse force. First of all, the authors usually 
consider point vortices. However, for small distances to 
the vortex axis, the corresponding velocity field cannot 
be taken literally since then the fluid velocity would di- 
verge and the Mach number of the problem would be- 
come infinite. It is not a priori clear that replacing the 
point vortex by a realistic vortex core would not affect 
the oscillatory motion of the vortex core due to the in- 
cident wave, and thus the momentum transfer between 
vortex and sound wave. Second, using simplified scalar 
wave equations instead of the full linearized hydrody- 
namic equations, one does not take into account the ef- 
fects due the interaction of the incident wave with the 
vortex core leading to its oscillatory motion accurately 
Q|. Finally, working, as is conventional, with an (un- 
physical) infinitely extended plane wave results in math- 
ematical difficulties (due to the non-convergence of sums 
over partial wave shifts in the scattering cross section 
E). We will show that no such singularities appear in 



the momentum transfer between sound wave and vortex, 
provided one is using a realistically shaped incident wave 
of a general form. 

In the following, we present a study of the scatter- 
ing problem which avoids using the simplifications dis- 
cussed in the above. In particular, we derive a global 
solution for the coupled hydrodynamic equations using 
the method of matched asymptotic expansions, matching 
the solution for the vortical region, with a realistic core 
of finite extension, smoothly to the far-field wave region 
solution, assuming a general incoming wave instead of a 
plane wave. As a result, we obtain a regular expression 
for the transverse momentum exchange between sound 
field and vortex. When interpreted in terms of two-fluid 
hydrodynamics, this momentum exchange exactly corre- 
sponds to the Iordanskii force. 



II. BASIC HYDRODYNAMIC EQUATIONS 

Our starting point are the Euler and continuity equa- 
tions for an ideal fluid described in terms of density p, 
velocity v and barotropic pressure functional P = P(p): 



| + v.(H = o, 

dv 1 
— + (v ■ V)v = — VP. 

at p 



(la) 
(lb) 



In order to consider small perturbations {pi,v\) on a 
given background solution (po,v ), which we require to 
satisfy the hydrodynamic equations, we introduce the fol- 
lowing expansions of velocity and density: 



v = v +vi, 



p = po+ pi- 



rn 



Inserting these ansatze into and equating terms of 
first order on both sides results in the following equations: 



dpi 
dt 



V ■ (pavi) + piV -v = 0, 



dvi , _ N _ / c 2 

— + (wi • V)v Q = -V — pi 

dt \p 



(3a) 
(3b) 
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where the sound speed c is as usual defined by c 2 = 
dP/dp and d/dt = (d/dt + v -'V) denotes the convective 
derivative. We now switch from (pi,vi) to another set 
of variables (V'jJil y i a the following decomposition of the 
perturbations |15| 



Po dy> 

c 2 di ' 



(4) 



The velocity perturbation is decomposed into a part with 
zero vorticity, V-0, where the pressure potential ip is reg- 
ular, and a part comprising the vorticity in the sound 
wave, 

Our linearized hydrodynamic equations 10 may then 
be rewritten as 

dK^)^ v ' MW+a (5a) 



df 



(£ • V)v - x w = 0. 



(5b) 



The above coupled pair of equations has first been de- 
rived by Bergliaffa et.al. using Clebsch potentials [l^ . 
Eq. jSa) with the vortical part of the fluid velocity per- 
turbations, £, set to zero, was derived in [lfiL Il7j. It can 
be shown [15], that in the case where the frequency fl 
of the sound wave is much greater than the background 
vorticity u>o = V x vq, there exist solutions where £ is 
smaller by a factor ujq/ £2 than ip. This leaves us with only 
one equation for the determination of the scalar quantity 
tp when the quasiclassical condition ojq/ ' fl <gc 1 is fulfilled; 
see the analysis of the wave region in the section to follow. 

We study an axisymmetric vortex, which may be de- 
scribed by the following ansatz for the background so- 
lution, employing a cylindrical system of coordinates 
(r, 4>, z) in what follows: 



Po = Po(r), v = v (r)e 4 



(6) 



The radial velocity vq is linked to the vorticity ujq = u>q e z 
(which we consider to be a fixed quantity) by the relation 



1 



vo(r) = - r'u>o(r') dr' . 



(7) 



We require the vorticity for large r to decay faster than 
any power of r. This will later on allow us to use an 
approach based on the method of matched asymptotic 
expansions [IsL Hfl| . For the velocity this assumption im- 
plies that, for r — ► 00, it is given by v = j/r (here 
and below we ignore exponentially small terms), with 
7 = J r'u)o(r') dr' , which is the circulation divided by 
2ir. While the continuity equation llat is automatically 
satisfied by the ansatz the density po must be de- 
termined such that the Euler equation ijlbll be fullfilled, 
which then leads to the following equation: 



Po 



Po vl 



(8) 



where the prime denotes derivation with respect to r. We 
immediately see that density variations are of second or- 
der with respect to the Mach number J( which is defined 
as jtfi = % /coo, with ^ being some typical velocity in 
the vortex and c x = c(r — ► 00) the sound speed for the 
fluid at rest at infinity. We require the Mach number to 
be a small quantity (M <C 1), which allows for an ex- 
pansion of the solution of the hydrodynamic equations in 
j$ in the section to follow. 



III. MATCHED ASYMPTOTIC EXPANSIONS 

To carry out the matching procedure, we divide the 
plane into two regions: An inner, vortical region charac- 
terized by r ~ , in which the vorticity is concentrated, 
and an outer, wave region, where r is of the order of the 
typical wavelength A = 2 / nf2/c Qa of the incoming wave, 
which we take within the matched asymptotic expansion 
procedure to be A ~ After having developed suit- 

ably scaled non-dimensional equations for each region, 
the solutions in both region are expanded in the Mach 
number and matched where the regions overlap. 

In the vortical region we rescale the spatial variable 
r according to r — > r and the velocity v Q according to 
Vq —>■ 9/vq. The typical time scale of the problem is 
given by the inverse of the frequency fi^ 1 ~ /'W . Con- 
sequently, the appropriate scaling for the time variable t 
is t — > (/%f)t, and the frequency Q takes the scaled form 
('W /)fi. Thus we have |ct>o|//2 = 0(1), which means, in 
the light of the remarks in the previous section, that we 
should replace £ by £/. With these scalings, and using 
Eq. JHJ) our Eq- ll5ajl becomes: 



v-(vv + €) 



— ¥ (d r 4> + ir 

dt z r 



(9) 

while Eq. Il5bjl stays formally the same. 

In the wave region, the appropriately scaled spatial 
variable is given by R = Jtr. Furthermore, we define 
the velocity potential \P in the wave region as &(R) = 
ip(r). Since in the wave region the velocity is one order in 
the Mach number smaller than in the vortical region, we 
introduce the wave-region velocity Vq — *M~ x vq = ^ e^. 
Now, in the wave region, Eq. © takes on the following 
form: 



(V 2 - d t > 



'(Vq ■ V)d t & + 0(Jt A 



(10) 



Here and in what follows, we use the convention that the 
nabla operator acting on a wave-region quantity corre- 
sponds to a derivative with respect to R. The equation 
for £ is not needed because we restrict ourselves to so- 
lutions where in the wave region the vortical part of the 
velocity perturbation, £, is exponentially small. 

In order to decouple the equations in the vortical re- 
gion, we split £ into an irrotational and a divergence free 
part according to 



€ = V/; - 



(11) 
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where r\ and £ are uniquely specified by demanding that 
77, — > as r — > 00. Our Eq. © then becomes 

V 2 (V> + »?) = (-ifl + v Q • V) 2 V 

- Jt 2 ^[d r (^> + 77) - r"%C] + 0(^ 4 )- (12) 

This Poisson equation may be solved with standard 
methods. Using a relation following from Eq. Il5bjl 



dt 



v x (V x £) + V(v ■ £) - (Vip + £) x w = 0, 



(13) 



Bessel functions for small i?, yielding 



.4 



l e i(t<t,-Qt) 



=±1 



V=±2 



-i(2t 



(18) 



This will be needed later to carry out the matching. 

In the vortical region, Eq. 111211 to this order gives the 
ordinary Poisson equation V 2 (ip^ + T)' )) — 0, with the 
following general solution 



and operating with the curl, we arrive at the following 
equation: 



(o) 



.,(0) 



■\ e i(e<t,-nt)_ 



(19) 



r dtp) ^ r dip 
d 

= -iv' Q — {ij + r 1 )-u Jo V 2 {tP + T]). (14) 

Note that here the fact that V • Vq = has been used. 
The above equation allows us to calculate £ after ip + 77 
has been calculated using ltl2|l . 

We now formally expand the general solution in the 
Mach number j# and then solve the resulting equations 
to the relevant order in the Mach number. In the vortical 
region the expansion takes on the following form 

ip = V-> (0) + .^ip (l) + Jf 2 i> {2) + ..., 

77 = ^°) +^V 2) +•••, (15) 

C = C (0) +^C (1) +^ 2 C (2) + ---; 



while in the wave region we have 

W = W (0) + + J( 2 >$ (2) + 



(16) 



We restrict our analysis to order 0{^ 2 ), because at 
higher orders the calculations become very involved. 



A. The solution to O(^# ) 

At this order, the solution in the wave region is deter- 
mined by the ordinary wave equation (V 2 + fc 2 )^ ) = 0, 
where the nondimensional wave number k is defined as 
k = (7. We expand the general solution to zeroth order, 
corresponding to our incoming wave, in the following way 
into partial waves of order £: 



<f (0) = J2 A,,J w {kR)e l 



(e<l>-nt) 



(17) 



In the wave region, r' e ' = 0(^f~' i '), ruling out the exis- 
tence of terms with \£\ > 0, since these would be of neg- 
ative order in the Mach number. Thus, ip^ + r)(°> = clq. 
A decomposition of into azimuthal modes accord- 



ing to C (0) = IX-00 d 0) M e i &'- at ') reduces Eq. <d to 



dr 2 



1 d 

r dr 



rfl — £vq 



C, (0) (r)=0. (20) 



We restrict ourselves to the case where there are no crit- 
ical layers present, fl 7^ £vo(r)/r for all £,r [13 ■ In 
this case, the solutions of this radial Rayleigh equa- 
tion are known to behave for large r as Q°\ r ) ™> 
oil + /?£? , ~^'), immediately giving at — for all I. 
Furthermore, using identity ltl3)l and vf^ — V?/>^ + 
= as well as Vtfj^ = — Vrj^°\ we can derive the 



equation — ifh/°) 



d<t> 



const., from which it fol- 
const. Therefore, as 



lows that = Ce l{ar / V0 ^ 
should be uniquely defined at each point in space, C=0 in 
the presently considered case; rf Q ^ thus equals a constant, 
which is zero because it must vanish at large r. 

Matching the solution in the vortical region with the 
wave region gives aQ — Aq, thus 



■iOt 



,(0) = £(0) _ 



0. 



(21) 



To this lowest (zeroth) order in the Mach number, there 
is thus no perturbative flow in the vortical region, but 
since ^ 0, there are density perturbations present. 



B. The solution to 



In the wave region, we have once more the ordinary 
wave equation (V 2 +fc 2 )<2^ 1 ) = 0. We require the solution 
to this equation to be causal; therefore we write 



where the Jut are Bessel functions. The form of in 
the vortical region may be determined by expanding the 



E 



B t Hf\kR) 



J(£<j>-Qt) 



(22) 
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where the Hankel functions of the first kind corre- 
spond to outgoing cylindrical waves. Due to the irregular 
behavior of the Hankel functions at the origin, none of 
the terms in the above sum matches to the solution in 
the vortical region, as becomes clear from Eq. J 2 fill below; 
therefore, B t = 0. Thus, <f (1) = in the wave region, re- 
flecting the simple fact that there is no sound radiation 
to linear order in the Mach number into the wave region. 

In the vortical region, + rf 1 ' is once more deter- 
mined by the ordinary Poisson equation, giving tjj^ + 
7]^ = J2e = _i be (At) 1*1 e l W- fi *) ; where the sum is lim- 
ited by the fact that higher-order terms would be of neg- 
ative order in the Mach number. The quantity ist 
determined by the following equation: 



-in 



v d 



V 2 



Using the ansatz = f(r) E*=±i he^-^ and 

V 2 «o = d r (r d r (rVo)) = u>' , we immediately see that 
the equation above is solved once one chooses / = — iv$. 
The asymptotics of the solution thus obtained for r — > oo 
is given by 



d_ 

~d4> 



<•(!) 



i=±i 



b-nt) 



(23) 



• 7 
-i — 

r 



!=±1 



i(itp-nt) 



(24) 



One gets the full solution to by adding the general 
solution of the homogeneous Rayleigh equation to 

(p \ Since Cp 1 ' —> for r — > oo, Q — ► for r — > oo as 
well. Using the homogeneous version of Eq. l|23jl . which 
can be transformed to a Rayleigh equation like Eq. (12011 ; 
the same argument as given after Eq. J20I) then yields 



Mi 



(1) 



0. Since we have 



form w ip e l ( k ' R ~ nt \ for which the expansion coef- 
ficients are given by Ae — ipoe' 1 ^ 1 * i' e ', with 4>k the an- 
gle which the propagation vector fc = fc(cos</>fc, sin <f>k,0) 
makes with the x axis. We therefore have e z x V^ 1 ' w 
e z x (^ip uJoek)e~' int . Then, using v « 5(^0 x r) near 
the origin, e z x (±ip uj e k )e- int w -^o(e" im e fc • V)v . 
We see that the term e z x VC'- 1 -' causes a displacement 
r (t) = ipo e ~~ inte k of the vortex line. It follows for the re- 
sulting displacement velocity JjT*o(i) = —iQipo e ~ int e k = 
V!f - ^ ^| T . = o, so that the vortex indeed moves with the ve- 
locity in the incident sound wave to this order in the 
Mach number. 



C. The solution to 0(JZ 2 ) 

In the vortical region, V 2 (ip^ + r/ 2 )) = —k 2 A Q e~ im , 
with the following solution: 



.2) 



(2) Me 



-iOt 



ce(kr)^e 



t=-2 



(27) 

Comparison with II18J1 leads to cq — c±\ — and c±2 — 
A±2/8. In the wave region, obeys a forced wave 
equation: 



(V 2 + k 2 )W^ = -2ikjR~ 2 d^°l 



(28) 



In the appendix, it is shown that the solution to this 
equation may be written as a contour integral. The result 
is (with C denoting the integration contour shown in the 
appendix): 



b-nt) 



2tt 



(29) 



where the expansion coefficients are given by 



e z x VC 



(1) _ il. 



]T b e (e lP + l le r y^-^ 



e=±i 



i(icj>-Qt) 



i=±l 



(25) 



for r — > 00, the asymptotics of rf^ for r — > 00 can be 
determined to be 77 W -► -2 ^=±1 eb t e l ^~ nt \ thus 
giving 



1> 



(1) 



60 



EM 



«=±i 



j(e<t>-nt) 



(26) 



for r — > 00. By comparison with the solution in the wave 
region, we obtain feo = and b±i = A±i/2. 

The rotational part of , e z x VC*- 1 -* , describes the os- 
cillatory motion of the vortex core induced by the sound 
wave [20], which may be seen as follows. Near the ori- 
gin, we may expand the incoming wave into a plane-wave 



MR) 



In t + i — ) exp 



kR f 1 

i^r \t - - 
2 V t 



dt 



(30) 



We can write the general solution to Eq. 128|l as ^ 2 ' = 



^p 2) + with being the general solution of the 

homogeneous equation. For large R, <Pp corresponds to 

(2) 

an outgoing wave (cf. the appendix). Therefore, tf^ 
should also contain only the outgoing wave: 

00 

^h 2) = E CeH { l\kR)e^- m \ (31) 

£=-00 

The remaining coefficients are determined by matching 
to the solution in the vortical region, giving Ci = for 
t ^ ±1 and 



,(2) 



>(2) 



C±i = ±ijk-fA ±1 . 



(32) 
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The solution to quadratic order in the Mach number, 
W^ 2 \ is thus made up of two parts. The first part, ^p 2 \ 
may be interpreted to represent the interaction of the 
incident wave with the long-range velocity field of the 

(2) 

vortex, while the second, W£ stands for the interaction 
of the wave with the vortex core [lflj . 



IV. FORCE ON THE VORTEX FROM PHONON 
SCATTERING 

In this section we will return to the original (dimen- 
sional) variables of section |HJ Far from the vortex, 
our solution takes on the following asymptotical form: 
ijj = Yli V^( r ) e 2 ^ -174 ' where, using the asymptotic form 
of the expansion coefficients &t (R) in the far-field solu- 
tion iJSnjl, derived in the appendix [Eq. l|A5Jl ]. we have 



e 2 



-At 



\j2itkr 

Here, the phase shifts 6e are given by 



{(1 + 2iS e )e ikr + e -i*r+iT\e\+if I _ 

(33) 



-sgu(*)§£ ;i?±l, 



(34) 



which are idential to those obtained by Sonin jj|. To eval- 
uate the momentum transfer between sound wave and 
vortex, we consider the time-averaged momentum-flux 
tensor, which reads 0,0, Q] : 

m \ ( c lo (Pi) ( v i)\ x , / v ton 

{ n ij) = I 2 Poc ~ ) ij + P°°\ VliVl ih ( 35 ) 

where p^ is the constant density far from the vortex. 
We have to insert into this expression the time-averaged 
expressions for the density and the velocity expressed in 
terms of ip and £ according to 10} . Since we are looking at 
the region far from the vortex, we may use the simplified 
expressions: 



Pi 



.^Re[^] = ^V-^), (36) 



2c c 



v u = Re [ditp] = -(di?p* + ditp), 



(37) 



which, when inserted into the time-averaged momentum- 
flux tensor Q35ll . yield 

(nij) = ^[Wdrf + diWrt* + (fc 2 |vf - |v^| 2 )%]. 

(38) 

This comprises a general expression for the time-averaged 
momentum-flux tensor in the far- region where vq ~ 0. 

The projection F^, of the force on the vortex per 
unit length F in the direction of the unit vector no = 
(cos(/>o,sin</> ,0) is given by 



F, = (F.n ) = -f(n ij )n i 6a j . 



(39) 



Inserting l|38j) into the above expression, and choosing a 
cylindrical surface around the vortex as our contour of 
integration, we get 



F, 



<t>0 



Po 



) r cos(</> — 4>o) (d r tp*d r t/j + k ip*^) 



(40) 

Making use of 13311 . and performing the angular integra- 
tion, we find 



= pook Re V W+1 - ]l+1] e^A}A l+1 {5 e - 5 t+1 ) 



(41) 

Employing the obtained phase shifts for sound-vortex 
scattering, Eq. l(3"4*j) . we arrive at the following final ex- 
pression for the force: 



-7 Re 



(42) 

The above equation represents the general result for the 
force on the vortex from phonon scattering, containing 
partial waves up to order \£\ = 2, valid for an arbitrary 
incident wave form. 

Note that there are no singularities in (jl^l) present due 
to partial wave summation, because we have assumed an 
arbitrary incoming wave form right from the start. Only 
now, after obtaining the general result for the force in 
terms of the expansion coefficients At, we use a plane 
wave approximation for the incoming wave near the ori- 
gin. Then, with A t = ijj e- te ' t ' k i^, 



Pook 2 



7r-y | Vo | 2 sin(0o - </>*). 



(43) 



Using the following expression for the momentum density 
of a plane sound wave, j — |?Ao| 2 ^—k Eq. 1(3311 may 
easily be generalized to yield the parallel and transverse 
components of the force, 



o, F ± = -rx j, 



(44) 



where r — 2irje z and j is the momentum density of 
the wave at the location of the vortex. The coordinate 
invariant form of iHH is obviously not restricted to the 
case where the vortex axis points into the z-direction. 
For a vortex in a superfluid we have to replace j by the 
thermally averaged momentum density of the phonons 
j p h, which in terms of the quantities of two-fluid hydro- 
dynamics reads j p h — Pn(v n ~ v s ), where v n and v s 
are normal and superflow velocities, respectively. For a 
vortex which is at rest with respect to the superfluid con- 
densate (vl = v s ) this gives 



-Flordanskii = pn[T X (v L - «„)], 



(45) 



which is the Iordanskii force. We remark that the argu- 
ments presented in [8] suggest that the above form of the 
force remains invariant (to the considered order in the 
Mach number) , also if the vortex moves with respect to 
the superflow. 
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CONCLUSION 



Eq. ffify reduces to 



We have solved the problem of scattering of sound 
waves by a vortex using an expansion into partial waves 
and the method of matched asymptotic expansions, up 
to the relevant quadratic order in the Mach number. To 
the best of our knowledge, the analytical form of the so- 
lution for the scattered wave in the form of a Schlafli 
contour integral, Eq. Ij29ll . has not been given yet in 
the literature. By evaluating the momentum-flux tensor 
far from the vortex, we derived a general expression for 
the force exerted by the sound wave on the vortex. The 
force depends only on the lowest partial waves (|£| < 2) 
and is thus completely determined by the shape of the 
incident wave in a small region around the vortex axis 
with a diameter of the order of a few wavelengths (an 
analogous result was obtained by Shelankov for an inci- 
dent wave of beam-shaped form (IS]). We demonstrated 
that the usage of a "physical" incident wave of arbitrary 
shape (and finite extent) removes the mathematical diffi- 
culties which plagued previous work employing the par- 
tial wave method, where the total transverse cross sec- 
tion, obtained after performing a summation over all par- 
tial waves, contained singularities (was not convergent) 
PJ- 0.01- The resulting global, singularity-free solution of 
the scattering problem we presented then leads, in con- 
junction with conventional two-fluid hydrodynamics, to 
the Iordanskii force on the vortex, confirming its exis- 
tence as a result of the asymmetrical scattering of sound 
waves at a vortex. 

Direct experimental confirmation of the Iordanskii 
force has so far been elusive. A conceivable experimen- 
tal procedure to verify the existence of the Iordanskii 
force in a superfluid is to study the influence of sound 
pulses in Bose-Einstein condensates, created by Bragg 
scattering 21], on the vortex. Using Bose-Einstein con- 
densates has, inter alia, the advantage of being able to 
create sound wave pulses in an accurately controlled man- 
ner, with variable large momenta of the order of c/7. We 
point out that our approach predicts the transverse force 
correctly also if the incoming wave is not a plane wave, 
which will in general be the case in the inhomogeneously 
trapped condensates. If the vortex is pinned by an opti- 
cally created external potential to fix it in one position, 
and, shortly after the sound wave has been created, the 
pinning potential is turned off, one should in principle be 
able to observe the transverse displacement of the vortex 
due to the sound wave scattered by it. 



APPENDIX A: FORCED WAVE EQUATION 



In this appendix, we derive the solution 129)1 to the 
forced wave equation Ip8f . Employing the ansatz (E 



W 



d 2 
dR 2 



R— + k 2 R 2 - £ 2 
dR 



<P e = -4\£\ni J w {kR). 

(Al) 

We write the solution to the above inhomogeneous Bessel 
equation as a contour integral. On the right-hand side, 
using the Schlafli integral representation for the Bessel 
functions [23, with the contour C shown in Fig. we 

express J\i\{kR) as 



J\i\{kR) = / <\\j> 



kR 
~2~ 



dt 



(A2) 



On the left-hand side, inspired by the above Schlafli in- 
tegral representation, we try the following ansatz: 



9i(t) exp 



kR 



dt 



(A3) 



After partially integrating the left-hand side and using 
the well known recursion relation Je-i(x) + Je+i(x) = 
■^rJi(x) for the Bessel functions on the right-hand side, 
we find that the function gg(t) is determined by = 
t~ l . If we choose the integration constant such that 



t (t) = hxt + i\ 



(A4) 



then the <Pi correspond to causal solutions of l)Al|l . This 
may easily be verified by applying the method of steepest 
descent, yielding the following asymptotics: 



>I>i 



-2tt 



(A5) 



(2) 

for R — > 00. Notice that ^ is regular at the origin, 

(2) 

% ' -> as R -> 0, as can easily be confirmed by ex- 
panding the integrand for small R and integrating. 



Im [f] ■ 




Re [tl 



FIG. 1: Schlafli contour of integration 
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